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0^ , Grover's quantum algorithm for an unstructured search problem and the Count 

I algorithm by Brassard et al. are generalized to the case when the initial state is 

Q"^ ■ arbitrarily and maximally entangled. This ansatz might be relevant with quantum 

subroutines, when the computational qubits and the environment are coupled, and 
in general when the control over the quantum system is partial. 



Abstract 



1 Introduction 



I In the recent years quite a significant progress has been made in the theory 

of quantum computation [1-3], both at the theoretical and experimental level 
[4-5]. In the quest for quantum algorithms, in particular, after the discovery by 
Shor [6] of an algorithm for factoring integers (which achieves an exponential 
speed up compared to the best classical algorithm currently known), one of 
the main successes has been Grover's algorithm for the unstructured database 
search [7]. Grover considered the problem to find a 'good' file, represented as 
the state \g >, out of files \a > ; a = . . . N — 1. The algorithm starts 
with the preparation of a fiat superposition of all states \a >, i.e.Q 

1 N-l 

\^,>=w\o>=^Y. I«>' (1) 



N a=0 



Assuming, without loss of generality, that is a power of two. 



Preprint submitted to Elsevier Preprint 



1 February 2008 



where W is the Walsh-Hadamard transform, and assumes that there is an 
oracle which evaluates the function H{a), s.t. H{g) = 1 for the 'good' state 
\g >, and H{b) = for the 'bad' states \b > (i.e., the remaining states in 
the set of all the a's). The unitary transformation for the 'search' of \g > is 
usually defined [7-8] in terms of the operator Gh = —WSqWSh, with the 
inversion operators Sq = I— 2|0 >< 0| and Sh = 1 — 2 J2g \9 >< S'l.Q Iterating 
Gh for n ~ 0[^/N] times on (1) then produces a state whose amplitude is 
peaked around the searched item \g >. Classically, it would take of the order 
of 0[N] steps on the average to find the same element g, so that Grover's 
quantum method achieves a square root speed up compared to its classical 
analogue. Subsequently, Grover's algorithm has been extended to the case 
when there are t 'good' items \g > to be searched or when the number of 
'good' items is not known in advance [8-11]. The number of steps required 
in these cases is of the order of 0[^N/t], again a square root improvement 
with respect to the classical algorithms. Then, Brassard et al. [12] combined 
Grover's operator and Shor's quantum Fourier transform in an algorithm that 
counts the number of 'good' items present in the flat superposition (1) with 
an exponential precision and a success probability exponentially close to one. 
It has then been shown that Grover's algorithm is optimal [13-14], that a 
single (complex) oracle query in Grover's algorithm might suffice in certain 
cases for finding the 'good' states [15-18], the algorithm has been exploited in 
A^P-structured search problems [19-20], its Hamiltonian formulation [21] and 
robustness discussed [22]. Although Grover's algorithm was originally devised 
assuming that the starting state is to be prepared in the flat superposition form 
(1), subsequent work [8, 11, 23] showed that this condition can actually be 
relaxed and that one can also work with the general initial pure state lipo >= 
U\0 >, i.e. replace W by an arbitrary unitary transformation U, and then 
use the operator Gh = —USoU~^Sh- It was also shown [24] that Gh can be 
interpreted in terms of a rotation of \ipQ > towards the vector representing the 
'good' states, \w >= J2g \g > /Vi, which is a product of two refiections in the 
two-dimensional space spanned by {ipo > and \w >. An explicit calculation 
for the case when the amplitudes of the initial superposition of states are 
arbitrary and unknown complex numbers was then made by Biham et al. [25], 
who found that one can still express the optimal measurement time and the 
maximal probability of success in a closed and exact form which depends 
only on the averages and the variances of the initial amplitude distribution of 
states.0 

^ In fact Sh can be easily implemented as an \a >-'controIled' unitary transforma- 
tion Uh hy tensoring \ipQ > with an extra ancilla qubit |e >= [|0 > — 11 >]/\/2, 
such that Unla > \e >= \a > |e + H{a) mod 2 >, thus obtaining Unltpo >— > 
-^^a=Q 1*^ ^ {—l)^^'^^\e >. Furthermore, the 'inversion about the average' oper- 
ator can also be compactly written as U^^ = WSqW = I — 2|-(/'o >< V'oI- 
^ Further discussion for the case of an arbitrary (non entangled) initial state and/or 
an arbitrary unitary transform in Grover's search can be found in Pati (1998) and 
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One of the main resources and ingredients of quantum computation lies, how- 
ever, not only in the possibility of dealing with arbitrary complex superposi- 
tions of qubits, but also in the massive exploitation of quantum entanglement 
(see, e.g., refs. [28-29]). In particular, it is not obvious to which extent one can 
directly use Grover's algorithm when the states to be searched are nontriv- 
ially coupled with the states of another system over which we cannot have a 
complete control, i.e. when the initial (normalized) superposition is given byQ 



where fa is an arbitrary mapping (not necessarily one to one). In fact, in 
Grover's algorithm the application of any unitary transformation acting on 
the computational states \a > would automatically affect also and nontrivially 
the \ fa >'s, by producing a complicated mixing of the original entanglements. 
An interesting case may arise, for instance, when the system \ fa > represents 
the environment. Of course, one might trace out (measure) the 'environment' 
and simply work with the \a > subsystem, which after the measurement of 
the \ fa > subsystem will collapse into a generic mixed state, i.e. an incoherent 
superposition of sectors of pure states (each characterized by a different value 
of fa) , and then simply study the dynamics in each of these sectors using the 
methods described, e.g., in ref. [25]. The problem with this approach is that 
sectors with different /a's will generally have different numbers of 'good' states, 
and thus require different computational time to find the 'good' states (in some 
cases failing, for the sectors where there are no 'good' states). One can give an 
upper bound on the success probability of the algorithm, which cannot exceed 
the expectation value of the operator laJ2g,g' \ fg >< fg'\, however, for generic 
situations in principle one could do much worse and there is no explicit scheme 
telling what one should expect. Alternatively, one might consider the setting 
in which Alice prepares the state (2), where the function fa is complicated 
and requires a lengthy computation, and sends it to Bob. Bob then wants to 
access the information contained in the \ fg >'s associated with certain marked 
states \g >, under the promise by Alice that there are t 'good' \g >'s, but 
without having an a priori 'local' knowledge about the /a's or how to compute 
them (except, possibly, the knowledge about certain 'global' properties such 
as the first order momenta of the \ fa >'s), and the simplest way to do that 
is to enhance the amplitude of the \g >'s using Grover's algorithm first and 
then read the associated \ fa >■ Finally, in number theory one might want to 
count (for example, when testing the Prime Number theorem, see ref. [30]) the 



Gingrich et al. (1999). 

^ We take the normalization conditions J2a=o < "I" >= < '^W ^aa', 
Ea=o < falfa >= N and Eg < fglfg ^^1, with A^i < iV and, in general, 

< falfb >+ Sab. 




(2) 
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number of primes k less than a given integer A*", and then need the quantum 
(entangled) superposition Ylk=o 1^ > >) where / is a 'primality flag' for 

k (i.e., / = 1 for k prime, and / = for k composite), which itself maybe the 
result of a previous lengthy quantum subroutine and which we want to reuse 
in further computations. 

The aim of the present paper is thus to generalize Grover's methods to the 
case when the initial state is given by the entangled superposition (2). 



2 Finding Good States 



The problem is to find the 'good' states \g >, promised to be in number 
t, from the initial normalized entangled superposition (2). By defining the 
remaining or 'bad' states as \b >, in number N — t, and applying Grover's 
unitary transformation Gh = —WSqWSh on the state > of eq. (2), it is 
easy to show (by induction) that the n-th iteration of Gh on > produces 
the state 



h\^>-4^{[T.\9> I/. >] + (-irE 1^ > \h >] 

VN I g ^ 

^|^>)|XW>+(^|6>)|F('^) >]j. 



2 
TV 



(3) 



Adopting a compact matrix notation, i.e. by substituting for \X^^^ >— > 
and >— Yn, writing ^„ = (X„,y„) and defining the matrix M = 

cos29{I + (Tx) + icTy (with the angle sin^ 6^ = t/N), the states > and 

ly^") > can be seen to satisfy the recurrence relations 



(4) 



which depend only on the number of 'good' items via the angle 9 and the 
following initial 'average' states 

1^(0) ^glfg > . 1^(0) >^ ^b\fb > 

' i ' ' N-t ^ ' 



via the quantities |C(") >= t\G^^^ > +{-lY{N - t)\B^^'^ >, with the sub- 
stitution IC*^"^-* >— Cn and writing C„ = C„(li 1) (subject to the the initial 
condition Xi = Yi = Ci). Eq. (4) can then be solved using standard tech- 
niques and, substituting back in eq. (3) we finally get for the n-th iteration of 
Gh on the entangled state \ip > 
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sin 2n6 ( 



sin 2nQ 



{-mh > 



(2sm^e\G^^^ > +(-l)"sin2e(tann^)(-i)"|S(°) > 



sin 2^ V 

^^[E\9>\ft^>+E\b>\t^ >]• 

V-'V g ft 



(6) 



As we can already note from expression (6), similarly to the case of the original 
Grover's algorithm acting on an initial flat superposition of states, also in the 
presence of entanglements Gff\ip > is periodic in n with period 7t/9, and one 
can anticipate that a Fourier analyis can still be performed in order to find an 
estimate of 9 (see next section) . 

It is then possible to show, following methods similar to those used in ref. [25] 

and after some elementary algebra, that the 'variances' of the distribution of 
the amplitudes of the initial entangled state are constants of the motion, i.e. 



2(n) , 



II d/i") > >) 



t 



= cr. 



2(0) 
G 



'G 



2{ — 



n) _ Eft II d/ft^"^ > >) f _ 2(0) _ o 



B-i 



N-t 

where || \x >|p=< x\x >, and jG*^""^ >= [Epug 
[Eft \fb^^ >]/ — t) are the 'averages' at time n. 



(7) 



>]/t and >= 



Using this fact, one can calculate the probability of picking up a 'good' item 
after n iterations of Gh over the initial entangled state >, defined by 
P(n) = E,</i'^)|/i"^>, as 



P{n) = Pav - APcos 2{2ne - 0^)6"^'^^ 
- Nal cos^ e 

< >+tan2^<G'(°)|G'(°) > 

where the complex angle = 0^ + i0/ is defined by the formula 

exp[2i0] 



Pav = 1 - AP - Nal cos^ 9 

AP = ^cos^^ 
2 



2<pf |Pi°)> 



[<Pf jPf > + <i^i°Vl°^ >1 



(8) 



(9) 



with |p|_ >= \B^^^ > ±itan^|G*^°^ >. As it can be easily seen from eqs. 
(8) and (9), the probability P(n) only depends on the first order momenta 
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of the distribution of the amphtudes of the initial entangled state, i.e. upon 
the initial 'averages' < G^^^IG^^^ >, < > and the initial 'variance' 

0"^, and on the number of 'good' items t. This is similar to the case of a non 
entangled initial state with arbitrary amplitudes [25]. The probability P{n) is 
maximized, Pmax = Pav + e^^'^' , at the times 

n, = [7T{2j + l)/2 + (j)n]/2e ; j e Z, (10) 

and the minimum number of iterations required to achieve the maximal prob- 
ability of success Pmax is uq. In particular, for the case of 6^ ~ \Jt/N <^ 1, 

we have that ~ 0{^N/t), which is the same as in Grover's original al- 
gorithm. Moreover, one can almost be certain to find a 'good' item \g > 
only provided that (for ^ -C 1) the initial 'variance' of the amplitudes of 
the 'bad' states is small enough, i.e., if Na'^ = e -C 1 then we have that 
Pmax — 1 — e ^ 0(1), independently of the values of the initial 'averages' 
< G^^^IG^^) > and < B^"^\B^"'^ >• Q A more detailed analysis of the pos- 
sible outcomes when the 'averages' and the 'variances' of the amplitudes of 
the initial entangled state \ip > are arbitrary and unknown (e.g., best average 
number of steps necessary to find the 'good' states \g >, etc..) can be done 
following the same lines as in ref. [25]. It is straightforward to show that the 
model of an initial state with arbitrary complex amplitudes described in ref. 
[25] can be recovered as a subcase of our algorithm provided one makes the 
substitutions \ fg >= ki\0 > and |/f, >= li\0 > (with ki and k arbitrary com- 
plex phases). This includes, in particular, Grover's original ansatz [7] for the 
choice ki = li = 1. 

Finally, it is also interesting to note that, for n = uj, we get ^/NG^lip >= 
Epk > [\fg > -(1 + (-l)^sin0^)|G(o) > +(-l)^cos0^cot0|fi(o) >]. Of 
course, unitarity prevents one to naively get only the exact contribution from 
the initial 'unperturbed' entangled states |/g > in |^ >• However, as some 
elementary algebra can show, it is still possible, for instance in the case of 
a large enough number of 'good' items g, i.e. for t/N < 0(1), to make the 
amplitude contribution coming from the other entangled states 1(5^°^ > and 



^ The probability of finding the 'good' states does not change with the number of 
trials n (besides the trivial cases of = it/2, t = N,0 = t = Oor, e.g., when the 
function /„ is one to one, since then the reduced density matrix pa = Trf\ip >< ^p\ = 
J2a 1^ X '^l) on which Grover's operator acts, is totally mixed, i.e. pa = I, and so 
invariant under any unitary transformation) when | < f!;^^\F^^ > | = 0, since then 
we get P{n) = Pav = 1 — ^ ~ const. The probabihty of finding the 'good' 

states can be small, i.e. P{n) ~ (5 ^ 1, if the initial variance o"^ ~ (1 — 6)/{N — t) 
((T^ ~ N~'^ for t <C N). Finally, the algorithm completely fails only when one of 
the previous conditions holds and there is a fine tuning Na'^ cos^ 9 = 1, for which 
P{n) = const = 0. 
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1^(0) > relatively small compared to that of \fg >, if j is even and provided 
that < (5(0)1(5(0) >~< 5(o)|B(o) >. 



3 Conclusions 



We have shown how to generahze one of the most useful algorithms discovered 
in quantum computing so far, i.e. Grover's algorithm for the search into an 
unstructured database, in the most general case when this algorithm is to 
be applied on a generic initial state in an unknown, arbitrary and maximally 
entangled superposition of qubits. This situation might arise, for example, when 
the computational qubits are nontrivially entangled with the environment, 
when one wants to search or count marked items in subroutines part of larger 
quantum networks, or more in general when part of the quantum system is 
not accessible. In particular, we have seen that even in this general case, 
the dynamics of the quantum entangled system is periodic, with the period 
depending on the number of 'good' items, and fixed by the first order momenta 
of the amphtudes of the initial state. Furthermore, the search algorithm still 
generically needs 0{^N/t) iterations to sort one of the 'good' items, and 
the maximum probability to obtain such 'good' items can be made close to 
one, provided that the initial 'variance' of the amplitude distribution of the 
'bad' states is small enough. In the appendix we also generalized the COUNT 
algorithm and showed that one can preserve a good success probability and 
a high accuracy in determining the number of 'good' items even if the initial 
state is entangled, provided that some conditions are satisfied by the 'averages' 
of the amplitude distribution of the initial state (for example, for the choice < 
(5(o)|(5(o) > > < 5(o)|5(o) > > 7rV(8^)). These results were not obvious 
a priori and constitute a non trivial feature for quantum computational states 
with entanglement. 
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A Counting Good States 



The algorithm COUNT, introduced by Brassard et al. [12] for the case of an 
initial flat superposition of states, essentially exploits Grover's unitary opera- 
tion Gh, already discussed in the previous section, and Shor's Fourier opera- 
tion F for extracting the periodicity of a quantum state, defined as^ 

1 k-l 

F|a ^^e^^^^/'^lo . (A.l) 



Then, the COUNT algorithm can be summarized by the following sequence 
of operations: 1) {W\Q >){W\Q >) = E™=o 1"^ > 'Ea=oW > 2) ^ 
(F ® /)[E™=o \m > G™(Ef=o' l« >)] 3) ^ measure \m >. 

The main idea at the core of the algorithm is that, since the amplitude of the 
set of the good states \g > after m iterations of Gh on |a > is a periodic 
function of m, the estimate of such a period by use of the Fourier analysis and 
the measurement of the ancilla qubit \m > will give information on the size t 
of this set, on which the period itself depends. The parameter P determines 
both the precision of the estimate t and the computational complexity of the 
COUNT algorithm (which requires P iterations of Gh [12]). 

Let us now discuss how one can use the previous results, obtained for the 
generalized Grover's algorithm when the initial state is arbitrarily entangled, 
to the case of the COUNT algorithm. We start by tensoring the state \ip > 
with logP ancilla qubitsQ set to |0 >, act on these qubits with W, obtaining 

>= Z]m=o 1"^ > 1"^ > /V^) with an \m >-'controlled' Grover operation 

on the state \ip > and with F on \m >, thus getting 

y-P-l 2mmn/P\ ^ 

1^/ ' G^\^ > . (A.2) 



As in the standard COUNT algorithm, requiring that the time needed to 
compute the repeated Grover operations is polynomial in logfc, leads to 
the choice P ~ O [poly {log k)] in eq. (A.2). Explicitly summing over n in eq. 
(A.2), after some elementary algebra we get 



0>\A> +\P/2 >\B>+^ 1^ > i<Pisi\C+ > 

m=0 



^ We assume, without loss of generality, that A; is a power of 2. 
^ Without loss of generality, we assume that P is a power of 2. 
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(A.3) 



where we have introduced the following quantities 



sin 7r(m ± /) 
Psin[7r(m±/)/P] 



m 



,i7r(m±/)(l-l/P) 



^, (A.4) 

7r 



with < f < P/2, and 1^4 >, |i? > and |C± > are certain mutually orthogonal 
states. 

At this point one can rewrite formula (A.3), similarly to how it is explained 
in ref. [12], in the general case when / is not an integer, distinguishing three 
possible cases. In particular, in the most general case in which l</<P/2 — 1 
we have 



\iPf > = !/-> \a^ > +\P -f-> \h > > |ci > 

+ IP-/+ > |di > +|Pi >, (A.5) 

where |Pi > is an 'error' term including all the other states in {ipf > not 
containing the ancilla qubits |/^ > and |P — > (with = [/] + 6f, 
= + ! and < 5/ < 1). After some easy algebra, one can show that 
the total probability amplitude in the first four terms is given by 

= [sm^e< G^^^IG^o) > +cos^e< p(°)|p(°) >]Si, (A.6) 



with Si = (■Sj+)^ + (-5^-)^ + (■Sp_j+)^ + (sp_^_)^, and it can be shown that 
S/tt^ < Si < 1. Similar calculations can be done for the cases when < / < 1, 
when the probability of obtaining any of the states |0>, |l>or|P — l>is 
given by W2 = {iVi + iV[(S2 - 1) sin^ 6 < G^^^G^^^ > +S2 cos^ 6 < pW|pW > 
]}/N (with S2 = (4)2 + (stf + (4_i)2 and S/vr^ < Si < 1), or when 
P/2 — 1 </< P/2, for which the probability to get any of the states |P/2 > 
o_r |P/2±1 > reads W3 = l-{A^i-iV[S3 sin^ 9 < G^jG^o) > +(S3-1) cos^ 9 < 
P(o)|P(o) >]}/N (with S3 = (4/2)2+(4/2-i)'+(4/2+i)'and8/7r2 < S3 < 1). 

The final step of the COUNT algorithm consists in measuring the first ancilla 
qubit in the state {i/jj >. As explained in ref. [12], if Wi > 1/2 (for i = 1,2 
or 3) then with high probability[^ one can still be able to find one of the 
ancilla qubits \f± > or |P - f± >, |0 >, |1 > or |P - 1 >, |P/2 ± 1 > or 
I P/2 >, respectively, for the three cases, and, therefore, evaluate the number t 

^ By repeating the whole algorithm many times and using the majority rule, Bras- 
sard et al. (1998), and eventually using R ancilla qubits |mi > ...|m/j > and acting 
with a '\mi > ...\mFt > -controlled' operation on the state \ip >, Carlini et al. 
(1999). 
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of 'good' states from sin 6 = \Jt/N and eq. (A. 4). For example, in the general 
case 1 < / < P/2 — 1, the condition that Wi > 1/2 is satisfied, e.g., for 
the choice of the initial 'averages' < G'-^^IG^^^ > > < > > 

7r^/(8v^)- Finally, although in general / is not an integer and the measured / 
will not match exactly the true value of / but give the approximate estimate 
i = A^sin^ ^(/)) the error over t for an entangled initial state will be |t — t| < 
7cN \n/P + 2Jt/N\ /P, i.e. the same as in ref. [12]. 
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